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A brief sketch of computer methods of involutivity analysis of differential equa- 
tions is presented in context of its application to study degenerate Lagrangian 
systems. We exemplify the approach by a detailed consideration of a finite- 
dimensional model, the so-called light-cone SU{2) Yang-Mills mechanics. All 
algorithms are realized in computer algebra system Maple. 



1 Introduction 

Among the properties of systems of analytical partial differential equations (PDEs) which 
can be studied without explicit integration there are two important ones: the question 
of their compatibility and the problem of posing of an initial-value problem (Cauchy 
problem), providing the existence and the uniqueness of an analytical solution. Both 
of these problems are crucial for the correct formulation of the evolution of degenerate 
Lagrangian dynamical systems. 

The main obstacle in the study of these problems for PDEs of a given order q is the 
existence of "hidden" integrability conditions. These conditions, q' < q order differential 
equations, are consequences of the given system of PDEs that can not be derived using only 
algebraic manipulations with the PDEs. The special class of PDEs, called an involutive 

*Talk given at the Seminar "Symmetries and Integrable Systems", Dubna, Russia, December 20, 2002. 



1 



system of PDEs, has all such integrability conditions incorporated in it. This means 
that differentiation of the system, called prolongation, do not reveal new integrability 
conditions. Examples of such involutive systems are the quasilinear systems of Cauchy- 
Kowalevskaya type (normal systems). 

The extension of a system by its integrability conditions is called a completion. From 
the completion point of view the linear homogeneous systems of PDEs with constant 
coefficients can be associated with systems of pure polynomial equations P0-[S1- The 
polynomial involutive systems ||4j provide a fruitful algorithmic tool in commutative alge- 
bra p. 

The general algorithmic foundation of the involutive approach is based on the concept 
of involutive monomial division invented in |3] and defined for a finite monomial set. Ev- 
ery particular division provides for each monomial in the set the self-consistent separation 
of variables into multiplicative and non-multiplicative ones. In the case of linear differ- 
ential system all its integrability conditions can be constructed by sequential performing 
of multiplicative reductions for non-multiplicative prolongations of the equations in the 
system 

The polynomial and linear differential involutive systems generate involutive bases of 
polynomial ideals and linear differential ideals jH], respectively. They are polynomial 
[7, and differential Grobner bases P HO] of special form. Though the involutive bases 
are generally redundant as the Grobner ones, their use make more accessible the structural 
information of the polynomial and differential ideals. The Janet bases j3] may be cited 
as typical representatives of involutive bases and have been used in algebraic and Lie 
symmetry analysis of differential equations [5]. 

The completion of differential equations to involution is the most universal algorithmic 
method for their algebraic analysis [5] and can be applied for the following purposes: 

• Check the compatibility of the systems of PDEs. In case of system inconsistency 
there is an integrability condition of the form 1 = which is revealed in the course 
of completion. 

• Detection of arbitrariness in the general analytic solution of analytic systems of 

PDEs HEiiini- 

• Elimination of a subset of dependent variables, that is, obtaining differential con- 
sequences of the given system, if they exist, which do not contain the dependent 
variables specified. 

• Posing of an initial value problem for a system of analytic PDEs providing existence 
and uniqueness of locally holomorphic solutions [21 E] • 

• Lie symmetry analysis of differential equations. Completion to involution of the de- 
termining equations for the Lie symmetry generators is the most general algorithmic 
method of their integration |^ I14j. 
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• Preprocessing for numerical integration. Certain non-commutative involutive bases 
can be used for automatic generation of finite-difference schemes PDEs |15| 

• Computation of the complete set of constraints for degenerated dynamical systems 
and their separation into first and second classes [T71 [TH| IT^ . 

Below the last application will be exemplified by studying the finite-dimensional de- 
generate system, the so-called light-cone Yang-Mills mechanics. In order to make presen- 
tation more transparent we start with a very short introduction of the main settings of 
the involutive method. 



2 Involutive polynomial bases 

The basic algorithmic ideas go back to M. Janet who invented the constructive ap- 
proach to study of PDEs in terms of the corresponding monomial sets based on the 
following association between derivatives and monomials: 



QfJ-l-l h/in^Q 

axf ■ ■ ■ dxf^r 



Note that the monomials associated to different dependent variables belong to different 
monomial sets. 

The association (^J) allows to reduce the problem of involutivity analysis of linear 
homogeneous systems of PDEs to the same problem for pure algebraic systems PP-|S],|S]- 
Having in mind this fact we shall state now the main definitions and results concerning 
the involutivity of algebraic systems. 

Let R = K[xi, . . . , Xn] be a ring of multivariate polynomials over a zero characteristic 
coefficient field K. Then a finite set F = {/i, . . . , f^} C -R of polynomials in i? is a basis 
of the ideal 

m 

<F>=<j\,...,U>={J2hifi \hjeR}. (2) 

i=l 

In the involutive approach to commutative (polynomial) algebra which is a map- 
ping of the involutivity analysis of linear PDEs |21 E] , for every polynomial in the finite 
set F the set variables separated into disjoint subsets of multiplicative and 

nonmultiplicative variables. 

To be self-consistent such a separation must satisfy some axioms j3] and every appro- 
priate separation generates an involutive monomial division in the following sense. Fix a 
linear admissible monomial order >- such that 

m 7^ 1 =^ m y 1 , (3) 
mi >- m2 '^=> mim >- (4) 

holds for any monomials (power products of the variables with integer exponents) m, mi, 
m2. Then for every polynomial f in F one can select its leading monomial lm{f) (with 
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respect to >~). All leading monomials in F form a finite monomial set U. If u & U divides 
a monomial w such that all the variables which occur in w/u are multiplicative for u, then 
u is called involutive divisor of w. We shall denote by L an involutive division, which 
specifies a set of multiplicative (resp. nonmultiplicative) variables for every monomial u 
in any given finite monomial set U and write u\lw if u is (L— )involutive divisor of w. In 
the latter case we shall also write w = uxv where, by the above definition, the monomial 
V = w/u contains only multiplicative variables. 

In the papers jH [TTl |T2] several involutive divisions were introduced and studied in 
detail. Here, as an example we present one called after M. Janet, who was one of the 
founders of the involutive approach to PDEs and who devised the related separation of 
variables p. 

Given a finite set U of monomials in {xi, . . . , Xn} and a monomial u = xf^ ■ ■ ■ x^" G U, 
a variable Xi {i > 1) is Janet multiplicative for u if its degree di in u is maximal among 
all the monomials in U having the same degrees in the variables xi, . . . , Xi-i. As for xi, 
it is Janet multiplicative for u if di takes the maximal value among the degrees in xi of 
monomials in U. If a variable is not Janet multiplicative for m in f/ it is considered as 
Janet nonmultiplicative. 

Consider, for example, a monomial set 

U = {xiX2, X2X3, xl} . (5) 

This gives the following Janet multiplicative and nonmultiplicative variables for monomi- 
als in U: 



Monomial 


Variables 


Multiplicative 


Nonmultiplicative 


X1X2 


Xl, X2, X3 






X2,X3 


Xl 


4 


X3 


Xi,X2 



Given a finite polynomial set F, a noetherian involutive division L, for instance, 
Janet division, and an admissible monomial order >-, one can algorithmically construct jl] 
a minimal L~involutive basis or L— basis G G Rof the ideal < F >=< G > such that for 
any polynomial / in the ideal there is a polynomial g in G satisfying lm{g)\Llm{f), and 
every polynomial g in G does not contain monomials having involutive divisors among 
the leading monomials of other polynomials in G. 

If F = {/i, . . . , fm} C -R is a polynomial set, L is an involutive division and >- is 
an admissible monomial order, then any polynomial p in R can be rewritten (reduced) 
modulo the ideal < F > as 



m 

p = h -^^CLijfi X Uij , 



(6) 



where elements (coefficients) of the base field K, L— multiplicative mono- 

mials for lm{fi) such that lm{f)uij ■< lm{p) for all i,j, and there are no monomials 
occurring in h which have L— involutive divisors among {lm{fi),...,lm{fm)}- In this 
case h is said to be in the L— normal form modulo F and written a.s h = NFl{p, F). 

If G is L— basis, then NFl{p,G) is uniquely defined ^ for any polynomial p. In this 
case NFl{p, G) = if and only if p belongs to the ideal < G > generated by G. Moreover, 
if the ideal is radical for which any its element (polynomial) vanishes at the common roots 
of all the polynomials in G if and only if this polynomial belongs to the ideal, then it 
follows that the condition NFl{p,G) = is necessary and sufficient for vanishing p on 
those common roots. 

It is important to emphasize that any involutive basis is a Grobner basis, generally 
redundant, and can be used in the same manner as the reduced Grobner basis [71IH]. 

The above described and some other properties of the involutive bases as well as the 
Grobner bases allow one to work fully algorithmically JHl with constraints in the 
case of degenerated dynamical systems of polynomial type. In particular, one can work 
algorithmically on the constraint surface. In the next section considering the light-cone 
Yang-Mills mechanics we shall show how the above mentioned ideas of the involutive 
analysis can be realized in the computer algebra system Maple. 

3 Application to Yang-Mills light-cone mechanics 

Before demonstration of computer calculations we formulate the Yang-Mills light-cone 
mechanics. At first we use the standard Dirac-Hamilton formalism for systems with 
constraints. Then we explain some computational aspects of deriving the same results 
implementing the Dirac's method in terms of involutive polynomial bases JT^ based on 
the Maple package. 

3.1 Dirac's constrained dynamics 

The Yang-Mills mechanics was formulated twenty years ago as instant form Yang-Mills 
field theory with spatially constant gauge fields and has been intensively studied during 
the last decades from different standpoints (see e.g [201-1201 and references therein). The 
light-cone version of Yang-Mills mechanics is formulated analogously, it follows from the 
classical action functional for Yang-Mills field theory with an Ansatz that the gauge 
potential is a function of the light-cone time only. So, we start with the action for the 
pure Yang-Mills gauge field in four-dimensional Minkowski space M4, endowed with a 
metric f] ^ 




I := — trF A *F, 

9 J Ma 



(7) 



^For other properties of the involutive bases, proofs and iUustrating examples see 
^In this paper we follow the notations of Ref. |27[ I28j 
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where is a coupling constant and the su{2) algebra valued curvature two-form 

F:=dA + AAA (8) 

is constructed from the connection one-form A. The connection and curvature, as Lie 
algebra valued quantities, are expressed in terms of the antihermitian su{2) algebra basis 
r° = a°'/2i with the Pauli matrices a" , a = 1, 2, 3, 

A = A V , F = F^t" . (9) 

The metric t] enters the action through the dual field strength tensor defined in accordance 
to the Hodge star operation = ^ y/r] e^ua/s F"^ ■ The light-cone coordinates = 

(x"*", x~,x^) are chosen as 



X 



± 



1 



(a;°±x^) , := (a;^A; = 1,2) (10) 

and the non-zero components of the metric t] in the light-cone basis are 

V+- = V-+ = -Vu = -V22 = 1 • (11) 
The connection one-form in the light-cone formulation is given as 

A = A+ dx~^ + A^ dx~ + Ak dx'' . (12) 

By definition the Lagrangian of the light-cone Yang-Mills mechanics follows from the 
corresponding Lagrangian of Yang-Mills theory if one supposes that connection one-form 
A depends on the light-cone "time variable" x~^ alone 

A = A{x^). (13) 

Using the definitions (|7j) and ()12|) we find the Lagrangian of the Yang-Mills light-cone 
mechanics 

L := ^ (F^ + 2 Fl, F\ - F^, F,%) , (14) 

where the light-cone components of the field-strength tensor are given by 

dA"" 

Fl_ = + e"*'^ A\ Al , (15) 
(9/1" 

Fl, = ^ + e'^''^A'^Al, (16) 

F-^ = ^-bc ^b_^c^^ (17) 
F-=,abc^b^c^ z,j,fc = l,2. (18) 

Choosing the light-cone coordinate x^ as an evolution parameter r 

r:=x+ (19) 
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we define the light-cone version of the Hamiltonian description of a system with La 

^aAdA 



grangian (fT^ . The corresponding Hessian is degenerate, namely corank||^^|| = 6, and 



thus the Legendre transformation ^ 

''a=W77 = 0, (20) 
7r- = -:^ = - (A"_ + e'^^^<AM , (21) 



dL 






= 0, 


dA% 


dL 


1 






dL 


1 


dAl 


^7' 



< = ^ = — e'^"" A" Al , (22) 
" dAl - 

impose the six primary constraints on the canonical coordinates. Thus the generalized 
Hamiltonian system is characterized besides the canonical Poincare-Cartan one-form 

Gc := vr+ dA% + tt" dA''_ + < dA'l - He dr (23) 

by set of primary constraints 

^«:=7r+ = 0, (24) 

Xr=/< + 6"'^AiA^, = 0, (25) 

satisfying the following Poisson brackets relations 

{^i'^^^'^} = 0, (26) 
{¥'i'^xU = 0, (27) 
{x-,x5} = -2^7V'^A^_r7.,. (28) 

In fl23|) the canonical light-front Hamiltonian is 

= y vt-tt; - e»^^ vr" + Al vr,^) + V(A,) , (29) 

where the potential term V is 

V{A,) = ^ [{AlA\) (A^A^,) - (A^A^,)] . (30) 

Following the Dirac formulation, the dynamics of the light-cone Yang-Mills mechanics 
is governed by the total Hamiltonian 

Ht = He- 2 ti (f/(r) if^'^) - 2 tr (Vfe(r) Xk) , (31) 



•^To simplify the formulas we shall use overdot to denote derivative of functions with respect to light- 
cone time r. Further, we shall treat in equal footing the up and down isotopic indexes denoted with 
a, b, c, d. 
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where U{t) and Vfc(r) are arbitrary SU{2) valued functions of the hght-cone time r. Using 

(2) 

this Hamiltonian we find that there are three secondary constraints (pa 

(^i2):=e,,,(A^_7r- + 4vre')=0, (32) 

obeying the so(3, R) algebra 

{^i^V^f^}=ea.c^i^). (33) 

Checking analogously the time evolution of the primary constraints xt have 

= xt = {xt, He} - 2 e^'^ At . (34) 

The analysis of this equation depends on the properties of the matrix Cab = e""^"^ At- 
This matrix is degenerate with a rank vary from to 2 depending on the point of the 
configuration space. If its rank is 2 then among the six primary constraints xt there are 
two first class constraints and maximum four Lagrange multipliers V can be determined 
from ()34|) . When the rank of the matrix Cab is minimal, the locus points are At = 
and all six constraints xt Abelian ones. For such an exceptional configuration the 
constrained system reduces to dynamically trivial one and hereinafter we shall consider 
the subspace of configuration space where rank||C|| = 2. For such configurations we are 
able to introduce the unit vector 

A"' 

N" = , , (35) 

^{Aiy + {Aiy + {Aty 

that is a null vector of the matrix || ^"-^'^ At ||, and to rewrite the six primary constraints 
xt as 

Xt^:=xt-{N\1) N\ (36) 
i^k:=N^xt- (37) 

Constraints xt± functionally dependent due to the conditions 

iV"xL = (38) 

and choosing among them any four independent constraints, we are able to determine 
four Lagrange multipliers V'lj_. The two constraints ipk satisfy the Abelian algebra 

{^,,^,} = 0. (39) 

(2) 

The Poisson brackets of the constraints tpk and (pa with the total Hamiltonian vanish 
after projection on the constraint surface (CS) 

{iJk,HT}\cs = 0, (40) 
{^^^\ Ht} \cs = (41) 
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and thus there are no ternary constraints. 

Summarizing, we arrive at the set of constraints (p^aK'ipk, ^^a \ Xk±- The Poisson bracket 
algebra of the three first ones is 

{^i'^^i'^} = 0, (42) 
{^,,^,} = 0, (43) 

{^i^^f} = ^a.cV^?), (44) 

{V^i^^ , i^k} = , ^} = , V^i^)} = . (45) 

The constraints Xk± satisfy the relations 

{X■^,X',A = -^9'^'""'A'LV^J, (46) 
and the Poisson brackets between these two sets of constraints are 

{V^a\xl±} = e'''^xl^, (47) 
W^:KxL} = {A.x',A = 0. (48) 

From these relations we conclude that we have 8 first-class constraints ip^a \ 'ipk, '^a'^ and 
4 second-class constraints xIl- This means that now constraints reduce the 24 constrained 
phase space degrees of freedom to 24 — 4 — 2(5 + 3) = 4 unconstrained degrees of freedom, 
in contrast to the instant form of the Yang-Mills mechanics where the number of the 
unconstrained canonical variables is 12. 



3.2 Computational aspects of Dirac-Grobner algorithm 

Now we shall discuss what kind of computer algebra manipulations are necessary to per- 
form in order to obtain the above stated results. We shall follow the general algorithm jTH| 
adapted to the computer algebra manipulations in theories with polynomial Lagrangians. 
This algorithm, called Dirac-Grobner algorithm, combines the constructive ideas of Dirac 
formalism for constrained systems with the Grobner bases techniques. 

Denote by qa and g^, 1 < a < 12, respectively, the generalized Lagrangian coordinates 
in dUl) listed as 

4^ A^ 4^ A^ A^ 4^ A^ A^ 4^ 4^ (AQ\ 

n.^, /i^, /i]^, /i^, /I2, /±2, ^-5 ^-5 ^- y^'-') 

and their velocities (time derivatives). Then the momenta are pa = < a < 12. 

To compute the primary constraints it suffices to eliminate the velocities qa from these 
system treated as polynomial in qa,qa,Pa- The elimination is performed by computing 
the Grobner basis [3 |H] for the generating polynomial set 

{Pa-^ I l<a<12} (50) 
oqa 

for an ordering (in Maplelexdeg) eliminating velocities qa- In the obtained set all alge- 
braically dependent constraints jH] are ruled out. 
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The canonical Hamiltonian (|29p is determined as reduction of 

PaQa - L (51) 

modulo the Grobner basis computed. Then computation of the Poisson brackets between 
the constraints is straightforward. 

Next step is the construction of the secondary constraints (j3!2|) . To find them we reduce 
the Poisson brackets of the primary constraints with the total Hamiltonian modulo the set 
of primary constraints. Again the Grobner basis technique provides the right algorithmic 
tool for doing such computations and to obtain a complete set of twelve algebraically 
independent constraints. In order to compare with the constraints given in the previous 
section, we represent them as 

^ = { ^fc, xl XU, a,b = 1,2,3; k = 1,2, (52) 

Next, to separate the complete set of constraints into first and second classes we compute 
the 12 X 12 Poisson bracket matrix on the constraint surface 

M:= a,(3=l,...12, (53) 

where faifp £ ^- Since the rank||M|| = 4 the complete constraint set T can be sepa- 
rated to four second-class constraints and eight first-class ones. To select the first-class 
constraints it suffices to compute the basis A = {ai, . . . , ag} of the null space of the matrix 
||M|| and then construct the first-class constraints as 

12 

First class constraints = ^^(ai)^/^ , l<z<8. (54) 

Q = l 

To extract the second-class constraints from one constructs 8x12 matrix || {a.i)a \\ from 
the components of the vectors in A and find the basis B = {bi, . . . , of the null space 
of the constructed matrix. Then every vector b G i3 yields a second-class constraint 

12 

Second class constraints = ''^^(hi) a fa, 1 < < 4 . (55) 

a=l 

As a result, one can organize the eight first-class constraints as ip^a\ipk-,^'a \ whereas 
the four algebraically independent constraints xL xl ^-re second-class. 

Relations (ji^ - lj^ revealing the structure of the gauge group generated by the first 
class constraints can also be computed fully algorithmically. To do this we extended the 
Maple package [TH] with a general procedure to represent the Poisson brackets of any two 
first-class constraints /„ and fp as a linear combination of elements in the set of first-class 
constraints: 

{fajp} = c\,f,. (56) 

With that end in view and in order to cope the most general case we implemented the 
extended Grobner basis algorithm [7j. Given a set of polynomials Q = {qi, . . . ,qm} 
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generating the polynomial ideal < Q >, this algorithm outputs the explicit representation 

ga = Kp qp (57) 

of elements in the Grobner basis G = {gi . . . , gn} of this ideal in terms of the polynomials 
in Q. Having computed the Grobner basis G for the ideal generated by the first-class 
constraints and the corresponding polynomial coefficients for the elements in G as 
given in (j^Tj), the local group coefficients c'^p (which may depend on the generalized 
coordinates and momenta) in (jSUj) are easily computed by reduction [71 |H] of the Poisson 
brackets modulo Grobner basis expressed in terms of the first-class constraints. 

However, the use of this universal approach may be very expensive from the com- 
putational point of view. For this reason our Maple package tries first to apply the 
multivariate polynomial division algorithm jS] modulo the set of first-class constraints. 
Due to the special structure of the primary first-class constraints that usually include 
those linear in momenta as in (j23), this algorithm often produces the right representa- 
tion ()57p: but unlike the extended Grobner basis algorithm does it very fast. Correctness 
of the output is easily verified by computing of the reminder. If the last vanishes, then 
the output of the division algorithm is correct. Otherwise the extended Grobner basis 
algorithm is applied. 

In our case the division algorithm just produces the correct formulas P^ - ()45|) for 
the Poisson brackets of the first-class constraints (p''a\ipk,'P^a^ ■ Similarly, one obtains the 
formulas (IIF)l)- (HH|l . 



4 The unconstrained system as conformal mechanics 

Now following the paper [27^ we demonstrate how using Hamiltonian reduction of the 
degrees of freedom one can derive the unconstrained version of the light-cone SU{2) 
Yang-Mills mechanics which coincides with the well-known model, the so-called conformal 
mechanics (2Mj. 

To show this equivalence we rewrite the model in terms of special coordinates as 
follows. At first we organize the configuration variables and A°i in 3 x 3 matrix Aab 
whose entries of the first two columns are A"} and the third column is composed by the 
elements A°i 

Aat:= (58) 

In order to find an explicit parameterization of the orbits and the slice structure with 
respect to the gauge symmetry action, it is convenient to use a polar decomposition for 
the matrix Aab 

A = OS, (59) 

where 5" is a positive definite 3x3 symmetric matrix, O(0i, 02,03) = e'^^'^^e'^'^'^^e'^^'^^ is 
an orthogonal matrix parameterized by the three Euler angles (0i,02,03) and 5*0(3, M) 
generators in adjoint representation {Ja)ij = ^iaj- Here we assume that the matrix A has a 
positive determinant and hence treat the polar decomposition (j^^ as a uniquely invertible 
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transformation from the configuration variables Aab to a new set of Lagrangian variables: 
six coordinates Sij and three coordinates 0j. ^ The polar decomposition induces 
the point canonical transformation from the coordinates Aab and Uab '■= ||7r"^ , tt"^ , 7r"~|| 
to new canonical pairs {Sab, Pab) and ((^a, Pa) with the following non-vanishing Poisson 
brackets 

{Sab , Pcd} = 2 i^ac ^bd + Sad ^bc) , (60) 
{(l>a,Pb} = Sab. (61) 

The expression of the old Uab as a function of the new coordinates is 

U = 0{P- kaJa) , (62) 

where 

ka = rab {Vb - ^bmn (SP)^,) , (63) 

7ifc = Sik — SiktrS and r]^ are three left-invariant vector fields on 50(3, M) group (the 
explicit form of 77^ in terms of the angular variables (pa , Pa can be found in |2ZI)- In terms 
of the new variables the constraints (pa and xf can be rewritten in the equivalent form 

Va=0, (64) 

Xai Pai ~l~ ^aij 1 jk ^kmni^SP^uin ~l~ ^amn Sm3 Sni (65) 

with vanishing Poisson brackets 

{Va,Xb^} = 0. (66) 

Passing to the new polar decomposition variables we achieve the complete separation 
of variables {Sab, Pcd), which are invariant under the gauge transformations generated by 
Gauss law constraints (pa , from the gauge variant one (0a, Pa). Owing to the separation 
to eliminate all gauge degrees of freedom related to this symmetry it is enough to project 
to the constraint shell described by the condition of nullity of Killing fields r/f'. The corre- 
sponding pure gauge degrees of freedom will automatically disappear from the projected 
Hamiltonian. 

In order to proceed further in elimination of remaining constraints we introduce the 
main-axes decomposition for the symmetric 3x3 matrix S 



S = R^{Xi,X2,X3)\ q2 \ Rixi,X2,X3), (67) 




^Note that the polar decomposition is vahd for an arbitrary matrix but the orthogonal matrix in (|59|l 
is uniquely determined only from the invertible matrix O = AS^^ , S = V AA^ and thus only in this 
case the polar decomposition (|59|l is a well-defined coordinate transformation. Concerning the degenerate 
matrices the more close and subtle analysis should be done. Here we note only that the set of n x n 
matrices with rank fc is a manifold with dimension k{2n — fc), but as distinct from the non-degenerate 
case now the manifold atlas contains with a necessity several charts. 
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with an orthogonal matrix R, parameterized by three Euler angles Xii X2, Xs- Because the 
Jacobian of this transformation is J (^^^^ j ~ Yll^b \ ^la ~ Qb \ the equation (jUTj) can be 
used as definition of new configuration variables: three "diagonal" variables (f/i, ^72, 93), 
eigenvalues of the matrix S, and three angular variables (xi, X2, Xs)) if only if all eigen- 
values of the matrix S are different. To have the uniqueness of the inverse transformation 
we assume here that qi < q2 < ■ ^ 

The momenta Pa and p^^, canonically conjugated to the diagonal and angular variables 
Qa and Xa, can be found using the condition of the canonical invariance of the symplectic 
one-form 

3 33 

Pab dSab = ^ Padqa + ^ Pxa ^Xa ■ (68) 
a, 6=1 a=l 0=1 

The original momenta Pab-, expressed in terms of the new canonical variables, read 
P = [PaC^a — - — - — C(a] R, (cyclic permutatious a ^ b c) . (69) 



a=l 



2 qb-Qc 



Here aa and aa denote the diagonal and off-diagonal basis elements for the space of 
symmetric matrices which obey the relations tr (aa^fe) = Sab, tr (aaOb) = 25ab, 
tr (aaab) = 0. 

The are three 5*0(3, M) right-invariant vector fields given in terms of the angles Xa 
and their conjugated momenta p^^ via 

41 = - sm xi cot X2 Pxi + cos Xi Px2 + Pxa ^ (70) 

smx2 

f, cosxi , , 

42 = COS xi cot X2 Pxi + sm Xi Px2 - ^xa > ( ' 1) 

^3"^ = Px. ■ (72) 
The constraints x rewritten in the main-axes variables take the form 



' ^ _ 1 1 ^ , 



T 

n 

a=l 



R, (73) 



where 



and Ur, = R. 



Pa = 1 ± - QaUaiqb ± ^c) , (74) 

qb ±qc 



a3 • 



^The variables qa in the main-axes transformation H67|l parameterize the orbits of the adjoint action of 
50(3, R) group in the space of 3 x 3 symmetric matrices. Whereas the consideration of the configuration 
Qi < 92 < 93 describing the so-caUed principle orbit class given below is correct, the treatment of all 
orbits with coinciding eigenvalues of the matrix S*, the singular orbits, requires more skillful treatment 
that is beyond the scope of the present paper. 
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As it was shown above the constraints Xi represent a mixed system of first and second 
class constraints, ipi and Xi± correspondingly. To perform the reduction to the constraint 
shell it is useful at first to introduce a gauge fixing condition and eliminate the two 
first class constraints ipi. The expression ()37p for the Abelian constraints ipi dictates the 
appropriate gauge fixing condition 

:=Ar-4» = 0, (75) 

which is canonical one in the sense that {tp^, tpj} = 6ij . The constraints ipi = together 
with the canonical gauge- fixing condition (j75p rewritten in terms of the main-axis variables 
fixes the canonical angular variables 

vr 

Xi = , Pxi = ' X2 = - , Px2 = . (76) 

Examine the remaining four second class constraints x} ^-^^d xl terms of the main- 
axes variables we find that the corresponding constraint shell can be described by the 
following conditions on the "diagonal" canonical pairs 

pi = 0, p3 = 0, (giigs)' = ±/ — • (77) 

Now using all above expressions for the constraints one can easily project the total 
Hamiltonian (jHT|) on the constraint shell and convince that the dynamics of the two 
unconstrained canonical pairs (^2,^2) and (xs^Pxa) governed by the following reduced 
Hamiltonian 

= T (^^ + i) • <^*^) 

where a = . Because the momentum p^^ is conserved one can identify the reduced 
Hamiltonian of SU (2) light-cone mechanics with the Hamiltonian of conformal mechanics 
whose "coupling constant" is determined by the value of a, while the gauge field coupling 
constant g controls the scale for the evolution parameter. 
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